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ABSTRACT 

We have measured i^s-band (2.0-2.3 /im) surface brightness fluctuations (SBFs) of 19 early-type 
galaxies in the Fornax cluster. Fornax is ideally suited both for calibrating SBFs as distance indicators 
and for using SBFs to probe the unresolved stellar content of early-type galaxies. Combining our results 
with published data for other nearby clusters, we calibrate /Cg-band SBFs using HST Cepheid cluster 
distances and /-band SBF distances to individual galaxies. With the latter, the resulting calibration is 

M Ks = (-5.84+ 0.04) + (3.6 ±0.8)[(V-I C ) - 1.15], 

valid for 1.05 < (V— I c )o < 1-25 and not including any systematic errors in the HST Cepheid distance 
scale. The fit accounts for the covariance between V — I c and Mk s when calibrated in this fashion. The 
intrinsic cosmic scatter of M k s appears to be larger than that of /-band SBFs. SO galaxies may follow 
a different relation though the data are inconclusive. The discovery of a correlation between Kg-band 
fluctuation magnitudes and colors with V — I c is a new clue into the star formation histories of early- type 
galaxies. This relation naturally accounts for galaxies previously claimed to have anomalously bright 
iGband SBFs, namely M32 and NGC 4489. Models indicate that the stellar populations dominating 
the SBF signal have a significant range in age; some scatter in metallicity may also be present. The 
youngest ages imply some galaxies have very luminous giant branches, akin to those in intermediate-age 
(few Gyr) Magellanic Cloud clusters. The inferred metallicities are roughly solar, though this depends on 
the choice of theoretical models. A few Fornax galaxies have unusually bright ifg-band SBFs, perhaps 
originating from a high metallicity burst of star formation in the last few Gyr. The increased spread 
and brightening of the /vg-band SBFs with bluer V— I c suggest that the lower mass cluster galaxies 
(< 0.1L*) may have had more extended and more heterogeneous star formation histories than those of 
the more massive galaxies. 
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1. INTRODUCTION 

Surface brightness fluctuations (SBFs) provide a pow- 
erful method to determine distances to early-type galax- 
ies (Tonry & Schneider 1988). Unlike many cosmological 
distance indicators, this method has a firm physical ba- 
sis: SBFs arise from Poisson fluctuations in the number of 
stars within a resolution element. This method uses the 
fact that the ratio of the second moment of the stellar lu- 
minosity function (LF) to the first moment of the stellar 
LF has units of luminosity. (See reviews by Jacoby et al. 
1992, Tonry 1997, and Blakeslee et al. 1999.) Using SBFs 
as a distance indicator requires that the bright end of the 
stellar LF in elliptical galaxies and spiral bulges is uni- 
versal, or that variations in the LF from galaxy to galaxy 
can be measured and calibrated. Tonry and collaborators 
have extensively pursued /-band SBF measurements (e.g. 
Tonry et al. 2000a,b). They have found /-band SBFs vary 
strongly between galaxies, but these variations are well 



correlated with V— I c galaxy color. 

Much less work has been done on SBF measurements 
in the near-infrared (1-2.5 /mi). Since SBFs in ellipticals 
are dominated by late-type giant stars, SBFs are brightest 
in the IR, with SBFs in /C-band (2.2 fj,m) about 4 mag 
brighter than in /-band (Luppino & Tonry 1993). The 
color of the night sky is even redder (J — K > 5) meaning 
IR observations need longer integrations to achieve the 
same S/N. However, the extreme redness of SBFs (V- 
K ss 7) compared to contaminating globular clusters (V— 
K « 2.5) means that SBFs can ultimately be measured 
in the IR to greater distances than in the optical, out to 
cz w 10, 000 km/s using 8-10 meter class telescopes (Liu 
& Graham 2001) or Hubble Space Telescope (HST; Jensen 
et al. 2001), and thus determine H at distances where pe- 
culiar velocities are negligible. IR measurements are also 
less subject to uncertainties in dust and extinction cor- 
rections than optical ones. Moreover, if IR SBFs can be 
shown to have a small intrinsic dispersion, they would be 
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a potent tool for peculiar velocity studies out to large dis- 
tances. 

As with any distance indicator, the utility of IR SBFs 
relics on the accuracy of their zero point and a thorough 
understanding of their precision. Published IR measure- 
ments to date (Luppino & Tonry 1993; Pahrc & Mould 
1994; Jensen et al. 1996, 1998, 1999; Mei et al. 2001) have 
all been in the K-h&nd. The current calibration of K- 
band SBFs rests on the work of Jensen et al. (1998). They 
establish if-SBF as a secondary distance indicator using 
observations of M31 and four Virgo cluster galaxies, with a 
Virgo distance from HST Ccpheid observations. They also 
use /-band SBF distances to derive a tertiary calibration 
for if-band SBFs from a larger set of 11 galaxies. 

Since SBFs depend on the second moment of the stel- 
lar luminosity function, they also provide data on unre- 
solved stellar populations unattainable from the first mo- 
ment alone, i.e. integrated galaxian light and spectra. For 
instance, IR SBFs are sensitive to the presence of inter- 
mediate age (t < 5 Gyr) asymptotic giant branch (AGB) 
stars. The SBF signal in the optical and IR is dominated 
by stars on the red (first ascent) giant branch and asymp- 
totic giant branch, /-band SBFs are predicted by models 
to be roughly degenerate to age and metallicity changes 
(Worthey 1993; Liu et al. 2000), hence their utility as a 
distance indicator and inadequacy for stellar population 
studies. On the other hand, IR SBF magnitudes (and 
optical/IR SBF colors) are predicted to show strong vari- 
ations with age and metallicity This can potentially of- 
fer insights into the star formation history of early-type 
galaxies, unique from ordinary integrated light and spec- 
tra. (We refer to Liu et al. 2000 for a full discussion of 
stellar population studies using SBF measurements.) 

In this paper we present results from a study aimed at 
(1) accurately calibrating IR SBFs for cosmological dis- 
tance measurements and (2) using multicolor SBF mea- 
surements in concert with state-of-the-art stellar popula- 
tion synthesis models to study the unresolved stellar pop- 
ulations of early-type galaxies. We have observed the K$- 
band SBFs of 19 early-type galaxies in the nearby Fornax 
cluster — Fornax is an ideal place to calibrate IR SBFs as 
it is the nearest cluster (d ~ 19 Mpc) that both is compact 
and has many ellipticals. The HST Key Project has mea- 
sured Cepheid distances to three galaxies in the cluster 
(Silbermann et al. 1999; Madore et al. 1999; Prosser et al. 
1999; Mould et al. 2000b). Fornax's small angular extent 
implies that errors due to cluster depth should be small, 
unlike the case for the comparably nearby Virgo cluster 
(e.g., Yasuda et al. 1997; Neilsen & Tsvetanov 2000). In 
addition, by observing galaxies with a spread in proper- 
ties (e.g., mean color), we can also measure any second 
parameter effects and the intrinsic IR SBF dispersion. 

In § 2 we describe our observations and basic reductions. 
Our methods for measuring the stellar SBF signal are pre- 
sented in § 3, and our results are presented in § 4. In § 5, 
we present the resulting calibration of the K g-band SBFs 
and discuss the stellar populations of the cluster galaxies, 
including evidence for recent star formation. For readers 
primarily interested in the results, we suggest skipping to 
§ 4 and 5. The results presented here supersede the pre- 
liminary versions published in Liu et al. (1999) and Liu 
et al. (2000). 



2. OBSERVATIONS 
2.1. Blanco ^-m: SBF Data 

We observed Fornax cluster galaxies using the Blanco 4- 
m telescope at Cerro Tololo Inter-American Observatory 
(CTIO) on 1997 November 08-11 UT and 1998 Novem- 
ber 11-13 UT. We used the facility near-IR camera CIRIM 
with the Ks-b&nd filter (2.0-2.3 /xm; McLeod et al. 1995). 
The camera employs a Rockwell NICMOS-3 256 x 256 
HgCdTe array with a pixel scale of 0'.'419 pixel" 1 (§ 2.4) 
when used with the f/8 secondary. Both runs were photo- 
metric, except for parts of the very first night and the end 
of the very last night. The typical seeing FWHM was Of! 8, 
with a full range of Of! 6 to l'.'l. 

Our sample is listed in Table 1. Integrations on each 
galaxy were interlaced with blank sky fields in pairs, with 
data taken in an ABBA pattern, starting with the sky field. 
Offsets from galaxy to sky fields were typically several ar- 
cminutes, with larger offsets for brighter galaxies. The sky 
fields were chosen to be free of bright stars and away from 
other large galaxies. For both the galaxy and sky posi- 
tions, the center of each field was ultimately dithered in 
a 2 x 2 square pattern over the detector which was then 
repeated. An off-axis CCD camera mounted on a motor- 
ized X/Y stage was used to guide the telescope during the 
integrations. When moving between galaxy and sky fields, 
the stage was moved to maintain the position of the guide 
star on the CCD camera so that the same star was used for 
both fields. Mindful of the evidence that low S/N can lead 
to large systematic errors in SBF measurements (Jensen 
et al. 1996), we made sure to acquire plenty of integration. 
The total integration on each galaxy was typically 32 min, 
with an equal amount of integration on the blank sky field. 

2.2. CTIO 1.5-m: Surface Photometry 

Because most of our targets fill the field of view of 
CIRIM on the Blanco 4-m telescope, we also obtained 
wider-field images using CIRIM on the CTIO 1.5-m tele- 
scope. CIRIM on the 1.5-m telescope with the f/8 sec- 
ondary has a plate scale of 1//14 pixel -1 with the Kg fil- 
ter. Data were obtained in a similar fashion as for the 
4-m data, using interlaced pairs of sky and galaxy fields. 
The 1.5-m observations were not guided. The galaxy was 
stepped in a 2 x 2 pattern over the detector to cover a final 
area of w 8' x 8' centered on each galaxy. Offsets between 
galaxy and sky fields were s=s 8' — 10'. The total integration 
on the center of each galaxy was typically 16 min. 

Our first 1.5-m run was 1997 October 11-12 UT. Con- 
ditions were photometric for the second half of the first 
night and most of the second night. Galaxies observed 
during non-photometric conditions on the first night were 
re-observed on the second night during photometric con- 
ditions. The atmospheric seeing was variable, ranging 
from 1'/ 5-270 FWHM. Our second 1.5-m run was 1998 
November 08-10 UT, and conditions were photometric for 
all the data presented herein. The seeing was variable, 
and the images have an angular resolution ranging of l'.'l— 
175 FWHM. 

2.3. Reductions 

The data were reduced in a mostly standard manner for 
near-IR images. We took great care to ensure that the re- 
duction process did not introduce any spatial correlations 
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on the scale of the point spread function (PSF), which 
would contaminate the SBF signal and lead to systematic 
errors. All reductions and analysis were done using cus- 
tom software written for Research System Incorporated's 
IDL software package. 

We constructed an average bias frame from the median 
average of many dark frames for every combination of in- 
tegration time and coadds. This average bias was then 
subtracted from all images with the same combination of 
integration time and coadds. We created a flat field for 
each night from a series of dithered images of the twilight 
sky. We used the twilight sky as it is expected to be the 
best approximation in the Ks-bemd to a surface of uniform 
illumination. However, since this bandpass also included 
thermal emission from objects at room temperature, the 
total illumination on the detector may not have been non- 
uniform due to, e.g., warm dust particles on the instru- 
ment entrance window or scattered thermal emission from 
the telescope structure. Therefore, a flat field created by 
direct averaging of the twilight sky images would have sys- 
tematic errors due to the non-uniform component of the 
illumination. Instead, we used an iterative-fitting scheme 
to solve for the relative quantum efficiency (QE) of each 
pixel using the entire set of images simultaneously. For 
each pixel (x, y) in frame i of a set of twilight images, we 
represented the measured counts DNi as 

DNi(x,y) = q(x,y)x Si + k(x,y) (1) 

where q(x,y) is the relative QE of the pixel (i.e., the flat 
field), and Si and k(x,y) are the number of counts from 
the uniform and non-uniform illumination components, re- 
spectively. By definition, the median values of q and k were 
1 and 0, respectively, and k(x, y) could be either positive 
or negative as it represented the deviations from uniform 
illumination. Also by definition, Si depended only on the 
frame number i, and since it came from the twilight sky, 
it was a monotonic function of i (decreasing for sunset 
twiflats and increasing for sunrise twiflats). Note that we 
assumed k(x,y) was time- independent, a reasonable as- 
sumption given that the twilight sky frames were taken 
over a period of only a few minutes. Equation 1 is lin- 
ear, so using the median count level as Si, we found the 
best-fitting q(x,y) and k(x,y). For each pixel, we used a 
cr-clipping scheme to flag and mask any outlier frames due 
to the presence of a cosmic ray or a star. The result was 
the flat field q(x, y) and a map of the non-uniform illu- 
mination k(x,y), which was of order a few percent of the 
mean twilight sky flux. 

After flat-fielding, the individual images still did not 
have a flat (uniform) background. This was corrected by 
sky subtraction, which served two purposes: (1) to remove 
(most of) the night sky flux, which can be >100x the mean 
galaxy surface brightness, and (2) to remove non-uniform 
structure arising from flat-fielding errors and the afore- 
mentioned thermal emission. Note that even if the flat 
fields were perfect, sky subtraction would still have been 
a necessary step because of (2). 

A preliminary sky subtraction of the images was per- 
formed to identify astronomical objects. Then for each 
image, we constructed a running sky frame from the aver- 
age of prior and subsequent images of the blank fields, ex- 
cluding any of the identified astronomical objects from the 
averaging. All the individual images, both of the galaxy 



and of the adjacent blank sky fields, were processed in this 
fashion. 

The sky emission changes temporally so the running sky 
frame needed to be accurately scaled before being sub- 
tracted. The usual way to scale the image is to multiply 
the local sky frame so that its median matches that of the 
object frame. However, we could not use this scaling ap- 
proach for the galaxy images — the scaling factor would be 
slightly too large since it would be based on the sum of the 
night sky brightness and the median galaxy surface bright- 
ness. The result would be a slight over-subtraction of the 
local sky frame. Instead, we scale the local sky frames to 
the median counts in unsubtracted galaxy images exclud- 
ing a large circular region centered on the galaxy, typically 
80-120 pixels in radius. For the largest galaxies, there was 
still be a very slight over-subtraction, but this error was 
corrected when we calibrated the DC sky level (§ 2.4). 
We used these same circular masks for the scaled sky sub- 
traction of the blank sky images. This ensured the blank 
sky images and galaxy images were reduced in basically 
an identical fashion, which was useful for the SBF power 
spectrum analysis (§ 3.3). 

We used the center of the galaxy to register the individ- 
ual frames, which were then shifted by integer pixel offsets 
and averaged to assemble a mosaic of the field. We did 
not employ any sub-pixel interpolation to align the images, 
since any such scheme would have introduced high spatial 
frequency correlations in the pixel-to-pixel noise and con- 
taminated the SBF measurements. Bad pixels were iden- 
tified by inter-comparing the registered individual images 
to find outlier pixels. These were then masked during the 
construction of the final mosaic. 

We observed the faint HST IR standards of Persson 
et al. (1998) as flux calibrators. Hence, our resulting 
magnitudes are Vega-based. On most nights, we observed 
standards over a range of airmasses to measure the zero 
point and extinction correction for each individual night. 
On some nights, it was not possible to observe standards 
over a wide enough range of airmass; in this case, we 
used the zero point from standards observed close in time 
and airmass to the galaxies and assumed an extinction of 
0.090 mag airmass -1 , typical for our other nights. Extinc- 
tion corrections were obtained from the the DIRBE/IRAS- 
derived values of Schlegel et al. (1998). 

2.4. Calibration of Surface Photometry 

The final step in the image processing was to photomet- 
rically calibrate the CTIO 4-m Ks-band galaxy images 
used for SBF measurements. We used the wider-field 1.5- 
m images to account for the DC sky background. First we 
extracted azimuthally averaged surface brightness profiles 
from the 4-m and 1.5-m images. In extracting the profile 
from the 1.5-m images, we used a 20" wide annulus lo- 
cated at the edge of the images, typically 3' in radius, to 
assess the sky level. We then compared the 1.5-m and 4- 
m profiles and simultaneously fit for both a multiplicative 
and additive offset between the two. The multiplicative 
term represented any zero point offset and the additive 
term was the sky background. In order to do this profile 
fitting accurately, we carefully calibrated the plate scale of 
the 1.5-m and 4-m data relative to each other, which was 
necessary for this profile matching scheme. (The absolute 
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two frames is 

x = (j/o — He) cos a + [xq — x c ) sin a 
V = (Vo - He) sin a - (x - x c ) cos a 



(3) 

(4) 



plate scale for the 1.5-m data was determined using stars in 
common with the Digitized Sky Survey.) The differences 
between the 4-m and 1.5-m photometric calibrations were 
quite small, typically < 1%, providing confidence in the 
independent calibrations. The fitted sky level was then 
subtracted from the 4-m images. 

We checked the accuracy of our photometric calibra- 
tions against the aperture photometry of Persson et al. 
(1979). For the six galaxies common to our two samples 
(NGC 1344, 1379, 1380, 1387, 1399, and 1404), our pho- 
tometry is on average 0.048 mag brighter than the Persson 
et al. large aperture (56" diameter) measurements, with a 
standard deviation of 0.025 mag. The Persson et al. data 
were taken with InSb single-channel photometer and the 
CIT if -band filter, which has a slightly redder bandpass 
(2.0-2.4 /iin) than our K$ filter. To estimate the expected 
offset in the photometry, we synthesized Ks~Kcit col- 
ors using IR spectra of solar-metallicity M giants (spec- 
tral type M0-M5 III) compiled by Pickles (1998), which 
should be representative of elliptical galaxy spectra for this 
purpose (Frogel et al. 1978). The synthesized Ks~Kcit 
color is —0.032 ±0.003 mag, which includes accounting for 
the filters' fc-corrections appropriate for the redshift of the 
Fornax galaxies. Therefore, we conclude the photometric 
calibrations of our images are in good agreement with the 
Persson et al. measurements. 4 

3. SBF MEASUREMENTS 

We measured SBF apparent magnitudes (m) in a simi- 
lar fashion to the method described by Tonry & Schneider 
(1988) and Tonry et al. (1990). The basic steps are: (1) 
measuring the galaxy's mean surface brightness profile, (2) 
cataloguing the globular clusters and background galax- 
ies to quantify their contribution to the fluctuations, (3) 
measuring the total fluctuation variance from the Fourier 
power spectrum of the model-subtracted image, and (4) 
determining m by subtracting the variance due to con- 
taminating sources from the total measured variance. 5 We 
now describe each of these steps in detail. 

3.1. Surface Brightness Fitting 

Elliptical galaxy fits were done using the isophotal fit- 
ting scheme from Williams & Schwarzschild (1979). Con- 
taminating point sources or galaxies were first flagged and 
excluded from fitting. Each isophote was represented with 
an ellipse 

x 2 + Sy 2 = A (2) 

where (x, y) are the coordinates for the ellipse in a refer- 
ence frame where the a;-axis coincides with the major axis 
and the isophote's center is at the origin. The quantity S 
is the square of the ratio of the major to the minor axis 
(1/5* = 1 — e 2 , where e is the ellipticity), and A is the 
square of the semi-major axis. If (xq, j/o) are the coordi- 
nates on the data frame, the transformation between the 

4 There is a larger discrepancy between our photometry and the smaller aperture (29" diameter) photometry of Persson et al. (1979). After 
accounting for the difference in filter bandpasscs along with the potential error from the centering of the photometer aperture and the relatively 
smaller offsets to blank sky positions, the remaining offset is f»0.1 mag. Using curves of growth measured from our images, the size of this 
offset suggests the 29" aperture of Persson et al. was in fact about 10% larger in diameter. Using a larger sample of galaxies, Pahre (1999) has 
also suggested the smaller apertures of Persson et al. were actually slightly larger than quoted. 

5 When referring to SBFs in this paper, we mean the fluctuations in the galaxy surface brightness which arise from the statistical variations 
in the stellar surface density, i.e., the stellar SBFs. Sources unrelated to the galaxy stellar population, unresolved globular clusters and faint 
background galaxies, also produce fluctuations in the galaxy surface brightness (see § 3.4). On occasion, we refer to the fluctuations and/or 
the variance due to these contaminating sources. But the use of "SBFs" implicitly means the stellar SBFs. 



where (x c , y c ) is the location of the isophote center in the 
data coordinates and a is the position angle of the ma- 
jor axis measured clockwise from the y^-axis. For a given 
isophote with a semi- major axis squared of Ai, we assumed 
the neighboring isophotes had basically the same shape. 
Thus, the only difference for different pixels close to the 
isophote was the flux F(x,y). We assumed the flux in the 
neighborhood of the isophote could be represented simply 
as 

F(x,y) = F t (A/Ay (5) 

where 7 is the local slope of the intensity distribution. 

At each semi-major axis, we solved for the flux Fi, axial 
ratio v£, position angle a, center of the ellipse (x c ,y c ), 
and radial dependence of the intensity profile 7. The 
fits were stepped outward logarithmically in radius, and 
all the parameters were allowed to vary as a function of 
semi-major axis. The fitted parameters were then spline- 
interpolated as a function of semi-major axis, and a galaxy 
model was constructed using bi-linear interpolation of the 
isophotal parameters. 

About half of the galaxies had obvious residuals from 
purely elliptical isophotes. In these cases, we also fit for 
the higher order harmonic content of the isophotes. The 
procedure we adopted is similar to the widely used scheme 
of Jedrzejewski (1987). At each semi-major axis, we first 
solved for the ellipse parameters, as described above. Then 
along this ellipse with flux Fi, we fitted the residuals as 

F(E) - F t = ^A n sm(nE) + ^ B n cos(nE) (6) 

n n 

where E is the eccentric anomaly and n > 3. Typically 
the n = 3 and n = 4 terms sufficed for a good fit, though 
some galaxies required up to n = 6. As is well-known for 
bright ellipticals, the amplitude of the cos(AE) term was 
the strongest of all the harmonics, which corresponds to 
boxy or disky isophotes. 

We note that in their SBF measurement procedure, 
Tonry et al. (1990) fit and remove large-scale residuals 
on spatial scales several times the size of the PSF after 
subtracting the galaxy model. The purpose of this step 
is to produce a very flat background, which is needed 
for identification of faint objects. These low spatial fre- 
quencies are therefore contaminated and excluded in their 
measurement of the SBF power spectra. We chose not 
do this step. The reduction of our IR data involved sky- 
subtraction, which is not the case for reduction of optical 
imaging data. Therefore, the background of our images 
is already very flat. When cataloguing the faint sources 
around the galaxy (§ 3.2), the SExtractor software did fit 
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the large-scale background in the images before identifying 
objects, but we did not use this fit for any other purpose. 
At any rate, we did not use such large spatial scales in 
fitting the galaxy power spectrum (§ 3.3). 

Figure 1 shows /Cs-band images for two of the galaxies 
in our sample. After subtracting the fitted galaxy model, 
the SBFs become apparent, along with the contaminating 
globular clusters and background galaxies. 

3.2. Measuring Globular Clusters and Background 
Galaxies 

3.2.1. Photometric Catalogs 

After subtracting the fitted model from the galaxy im- 
age, we compiled a photometric catalog of all astronomical 
objects using the SExtractor software of Bertin & Arnouts 
(1996), version 2.0.21. We multiplied each pixel by the 
square root of its exposure time to create a mosaic with 
uniform noise over the entire field. Objects in the noise- 
normalized mosaic were then identified as any set of con- 
tiguous pixels with a total S/N > 4 within an area equiv- 
alent to a FWHM-diameter circular aperture. Aperture 
photometry was then done on the original /Ts-band im- 
age. We used the resulting "MAG_BEST" magnitudes, which 
for uncrowded objects use apertures determined from the 
moments of each object's light distribution. This method 
is similar to that of Kron (1980) and is designed to recover 
most of an object's flux while keeping the errors low. 

We ran Monte Carlo simulations to quantify the com- 
pleteness of the resulting object catalogs. These were also 
done to verify the accuracy of the magnitudes measured 
by SExtractor. We inserted 30 artificial stars of known 
magnitude into the images and then processed the images 
with SExtractor in the same fashion as the original data. 
For the artificial stars, we used the same stars as those 
used for the SBF power spectra measurements, described 
below. For each magnitude bin, we ran the simulation 
30 times for a total of 900 stars. We compared the recov- 
ered magnitudes to the input ones in order to determine 
the random and systematic photometric errors. The arti- 
ficial stars were chosen to span a relevant range in magni- 
tudes in 0.25 mag steps. Typical 50% completeness limits 
for the galaxy images were K rj 19.5 mag. 

3.2.2. Measuring the Luminosity Functions 

Globular clusters and background galaxies fainter than 
the detection limit appear as unresolved point sources, 
which contribute to the fluctuations. We quantified their 
contribution using simple analytic approximations of their 
luminosity functions. Background galaxy counts are well 
described by a power-law distribution 

n gal (m) = A 10-"" (7) 

where n(m) is the number of galaxies per unit area 
per magnitude. We used the galaxy counts measured 
by Saracco et al. (1999) as these offer the best com- 
bination of depth (K w 17 — 22) and area coverage 
(20 arcmin 2 ). A weighted fit to their tabulated counts 
gives 7 = 0.40 with A = 10 473 deg -2 mag -1 for K = 21, 
or one galaxy arcmin -2 mag -1 at K — 18.05. The slope 



and amplitude of these counts agree well with deeper 
(K < 24) smaller-field counts of Bershady et al. (1998). 
The slope of the counts is somewhat larger than that from 
the deep imaging of Djorgovski et al. (1995), but the ampli- 
tude agrees reasonably well provided that the magnitude 
corrections advocated by Bershady et al. are used. 6 

The globular cluster luminosity function (GCLF) is well- 
described by a Gaussian function 



^glob 



(to) 



Na 



-(m-m°) 2 /2<T 2 



(8) 



where No is the total number of globulars, a is the width 
of GCLF in magnitudes, and m° is the peak (or turnover) 
magnitude. To accurately fit all three free parameters re- 
quires a large sample of globular clusters reaching to well 
past the turnover magnitude. Our if-band data contained 
relatively few globular clusters because of the bright sky 
background and the intrinsically blue color of metal-poor 
star clusters. The resulting if-band LFs reached only a 
few magnitudes above the peak of the GCLF so any fits for 
all three parameters would have been poorly constrained. 
Therefore, like previous optical and near-IR SBF studies, 
we chose values for m° and a and fitted only for the am- 
plitude of the GCLF. 

Four of the galaxies in our sample have high quality 
F-band data which reach fainter than the GCLF peak: 
NGC 1344, 1380, 1399, and 1404 (see compilation in Fer- 
rarese et al. 2000a). For these, we used the measured a and 
location of the GCLF peak, converting F-band m° into K- 
band by assuming V — K = 2.28, the average of the Milky 
Way and M31 globular clusters (Barmby et al. 2000). For 
the remaining galaxies, we assumed a = 1.4 mag, which is 
a good match for most giant ellipticals (Whitmore 1997; 
Blakcslcc et al. 1997; Harris 2000). Some of these galax- 
ies do have GCLF measurements, but they do not reach 
well past the GCLF turnover; hence, they were not well- 
suited for a determinations (e.g. Kohle et al. 1996). For 
the GCLF turnover magnitude, we used relative /-band 
SBF distances from Tonry et al. (2000b) to tie the galaxies 
with HST GCLF data to those without, incorporating the 
errors in the /-SBF distances into the LF fits. From the 
four Fornax galaxies with high-quality V-band data, we 
found a weighted average of < M^ >= —10. 04 ±0.10 mag. 
For the other galaxies, we used their /-band SBF distances 
to compute the expected m° K . Note that we were in fact 
using the relative /-band SBF distances to determine rrfi K ; 
thus, we are not directly dependent on the choice of the 
/-band SBF zero point (e.g., Tonry et al. 2000a). Also, 
as opposed to using/including Virgo cluster measurements 
(e.g. Seeker & Harris 1993), the use of Fornax galaxies 
alone circumvents concerns about differences in My be- 
tween Fornax and Virgo (Ferrarese et al. 2000b), perhaps 
due to environmental effects (Blakeslee & Tonry 1996). 



3.3. 



Fitting the Galaxy Power Spectra 

We first selected a region of the galaxy to analyze. 
We constructed a software mask defining an annular re- 
gion and excluding contaminating sources (globular clus- 
ters and background galaxies). The model-subtracted 

6 Our adopted parameters for the faint galaxy population differed from Jensen et al. (1998), who used galaxy counts from Cowie et al. (1994). 
It can be shown (from equation 9 of Blakeslee & Tonry 1995) that we consequently derive a residual galaxy variance 10° 1 ( mc "'~ 18.5) i ar g er) 
where m cu t is the cutoff magnitude. For m cu t = 19.5, this factor is 1.25. 
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galaxy image is then multiplied by this mask, and the 
Fourier power spectrum of the central 256x256 pixel re- 
gion was determined. For most galaxies, we used an an- 
nulus spanning 4"-24" in radius. For the faintest galax- 
ies (NGC 1336, 1373, 1375, 1380B, and 1419), we used a 
smaller annulus of 4"-16". Some of the galaxies had strong 
residuals in the their inner regions (e.g., small disks) which 
were not well-fitted by elliptical isophotes; for these, we 
chose a larger inner radius and increased the outer radius 
to maintain the same analysis area. 

At this stage, the image was composed of fluctuations 
due to the galaxy's stellar population (i.e., SBFs), fluctua- 
tions from the unresolved contaminating sources which are 
too faint to be catalogued, and Poisson shot noise from the 
sky background. Read noise was negligible since the IR sky 
background was very bright. In the absence of seeing, all 
the pixels in the image would be uncorrelatcd. Therefore, 
the power spectrum of the image would be white noise, 
i.e., equal power at all wavenumbcrs. The amplitude of 
the power spectrum would just be the total variance of 
fluctuations and sky noise. 

However, because the fluctuations are astronomical 
sources, they are convolved with the point spread func- 
tion (PSF) before reaching the detector; this introduces 
a small-scale spatial correlation. In contrast, the Poisson 
sky noise is uncorrelated since it is not convolved with the 
PSF. Therefore, the power spectrum of the galaxy image 
has two components: (1) rising power at low wavenumbcrs 
due to the PSF-convolved fluctuations and (2) white noise 
due to the shot noise. Notice that atmospheric seeing has 
allowed us to easily distinguish between these two compo- 
nents in the power spectrum. The galaxy power spectrum 
P(k) can then be represented as 

P(k) = P x E(k) + Pi, (9) 

where Pq is total variance due to the fluctuations, Pi is 
the white noise component, and E(k) is referred to as the 
expectation power spectrum. As we explain below, this is 
very nearly the power spectrum of the PSF. 

To measure the amplitude of Po, we first measured the 
PSF power spectrum. We used the star/galaxy classifi- 
cations produced by SExtractor to identify bright PSFs 
in the galaxy images. A box of 25" was used to extract 
the PSF stars, and other stars within this region were re- 
moved. We then normalized the PSF to unity flux. For 
9 out of 19 galaxies, multiple PSFs were available, and 
we used the weighted average of the resulting SBF magni- 
tudes for the final result. Comparing SBF measurements 
for the galaxies with multiple PSFs led to an error esti- 
mate of 0.08 mag from PSF uncertainties, which included 
errors in the PSF normalization. For galaxies with only a 
single PSF, we added this amount of error in quadrature. 

To properly model P(k), we needed to account for two 
additional factors aside from the PSF power spectrum: (1) 
the radial variation in the amplitude of the SBF variance, 
and (2) the effect of the software mask on P(k). First, in 
the galaxy image, the SBF variance depended on radius 
since the variance scales linearly with the galaxy surface 
brightness. To account for this we created a "window func- 
tion," which was the mask times the square-root of the 
galaxy model. The window function was proportional to 
the rms fluctuations from SBFs in the masked galaxy im- 
age. Next, we needed to properly account for the effect of 



masking the image on the galaxy power spectrum. It can 
be shown (Liu 2000) that the appropriate way to do this 
is to convolve the power spectrum of the window function 
(which contains the mask) with the power spectrum of the 
PSF. This convolution is then E{k), the expectation power 
spectrum. The net result is that the PSF power spectrum 
is slightly broadened by the mask. Once we had E{k), we 
simply fitted P(k) to solve for Po and Pi in Eqn. 9. 

The very lowest wavenumbers were not suitable for SBF 
analysis: flat-fielding and sky-subtraction errors produce 
extra power on the largest spatial scales which contami- 
nate the SBF signal. To assess which wavenumbers to ex- 
clude, we examined the power spectra of the reduced blank 
sky fields. As described in § 2.3, for each galaxy, observa- 
tions of the sky field images were interlaced with the galaxy 
images and used the same dither pattern. Also, the sky 
fields were reduced in a nearly identical fashion. The prin- 
cipal difference was that the frames used to construct each 
local sky frame for sky subtraction were on average farther 
separated in time than those used for sky-subtraction of 
the galaxy frames (this of course is inevitable) . We masked 
the point sources and galaxies in the reduced mosaics of 
the sky fields and then examined the power spectra. We 
found that the unusable wavenumbers were k < 16, cor- 
responding to spatial scales larger than about 1/16 of the 
detector size. Thus in fitting the galaxy power spectra, 
we used only wavenumbers from k = 20 to k — 128, the 
Nyquist frequency. 

Figure 2 shows the A"g-band fluctuation power spec- 
tra for our 19 Fornax galaxies and the resulting fits. The 
power spectra were two-dimensional images, and we did 
the fit with the full images. To show the results, we have 
plotted the azimuthally averaged one-dimensional power 
spectra. 

In order to gauge the errors in measuring Pq, we ana- 
lyzed the power spectrum of each quadrant of the galaxy 
image independently. We then used the average of these 
four fits to obtain Pq, using the number of unmasked pixels 
in each quadrant as the relative weights and the standard 
error of the fits as the uncertainty. We ran a series of 
Monte Carlo tests to determine if these quadrant-derived 
errors were reasonable. Using the surface brightness pro- 
file measured for target galaxies, we created 100 images 
of artificial galaxies, added SBFs with known amplitude, 
convolved with a Gaussian PSF, and added Poisson shot 
noise appropriate for our observations. We then fitted 
the resulting fluctuation power spectra. Over a range of 
S/N appropriate for our data, we found errors from this 
quadrant-averaging were accurate. 

For NGC 1389, the quality of the power spectrum fit 
appears to be somewhat poor. The S/N in the PSF used 
for the fit was good, so potentially the problem arises from 
a mismatched PSF. However, with only one PSF available 
in the field, we cannot directly determine if this was the 
origin of the problem. An independent observation of this 
galaxy would be useful. 

For NGC 1373 and NGC 1375, the PSFs used for the 
power spectra analysis were somewhat lower S/N than 
the rest of the sample. To gauge if the PSFs were suf- 
ficient to accurately measure m, we used the aforemen- 
tioned Monte Carlo tests to simulate the effect of noisy 
PSFs. We found for the PSF and galaxy power spectra 



Liu, Graham, & Chariot 



S/N of these two cases, the errors in the m results are 
expected to be negligible. For NGC 1380B, no suitable 
PSF was present in the galaxy images. We chose to use 
a bright star in the sky field for this galaxy as the PSF; 
this had a comparable FWHM to faint point sources in the 
galaxy image. (The above-described Monte Carlo simula- 
tions also demonstrated that systematic errors in m mea- 
surements scale roughly linearly with the mismatches in 
PSF FWHM.) 

3.4. Computing the Stellar SBF Variance 

The variance measured from the power spectrum is the 
sum of the variance from the galaxy's stars, which is the 
signal we desire, along with the variance from unresolved 
background galaxies and globular clusters. Having fitted 
the LF for these contaminating sources (§ 3.2.2) using an- 
alytic functions, we integrated these functions to compute 
the variance below the cutoff magnitude. 

For a power-law galaxy luminosity function, Blakeslee 
& Tonry (1995) show the variance per pixel due to the 
unmasked galaxies below the cutoff magnitude is 
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where p is the pixel scale in arcsec/pixel, m\ is the zero 
point of the image (magnitude of a star which produces 1 
DN for the integration time), m cut is the cutoff magnitude, 
and m g is the magnitude where the galaxy surface density 
is 1 mag -1 arcsec -1 . All the parameters in the equation 
were fixed except for m cu t, which we determined from the 
completeness experiments (§ 3.2). 

Similarly, having fitted for the normalization No of the 
GCLF, we computed the variance per pixel due to the 
globular clusters below the detection limit. Blakeslee & 
Tonry (1995) show that this is 
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where m\ is the zero point of the image, vrP K is the /T-band 
GCLF peak magnitude, a is the GCLF width in magni- 
tudes, and erfc(x) is the complementary error function. 
When calculating errors for <7q C , we included the formal 
fitting errors on No as well as the errors in the /-band SBF 
distances to the galaxies which governs the value of m° K 
used in the GCLF fitting. 

The fitted amplitude of the power spectrum, Pq (eqn. 9), 
was then corrected for the variance P r due to the unre- 
solved globular clusters and galaxies. The quantity P r was 
the sum of the galaxy and globular cluster variances (p 2 „ a i 

and o~q C , respectively) divided by the mean galaxy surface 
brightness per pixel in the fitting region. This contami- 
nating variance amounted to about 10-30% of the total 
variance. The remaining quantity (Pq — P r ) was then the 
variance due solely to the stellar SBFs. This was converted 
to an apparent magnitude, rn,K s > using the measured pho- 
tometric zero point and extinction correction. The final 
errors on rn~K s comprised the quadrature sum of the errors 
in the photometric calibration (0.02 mag), the PSF uncer- 
tainty for galaxies with only a single PSF (0.08 mag), and 
the measurement errors in P r and Pq. 



4. RESULTS 

Table 2 summarizes our measurements. For each galaxy, 
we tabulate the mean /Cg-band surface brightness in the 
measurement region, the ratio of (Po — P r )/Pi as a mea- 
sure of the S/N of our data, the final /Cg-band SBF appar- 
ent magnitudes (rn~K s ), the /iTg-band SBF absolute mag- 
nitudes (Mk s ) derived from /-band SBF distances to in- 
dividual galaxies and a Ccpheid distance to the Fornax 
cluster, and the (distance-independent) I c —Ks SBF color. 

We have five galaxies in common with the /C'-band (1.9— 
2.3 /jm) measurements of Jensen et al. (1998). The overall 
agreement between our m's is good, with an average dif- 
ference of 0.12±0.04 mag (0.07±0.06 mag if we exclude 
their low S/N data for NGC 1339), in the sense that our 
measurements tend to be fainter. Our data use the Ks- 
band (2.0-2.3 /um) filter, and the difference in filter is ex- 
pected to make the Jensen et al. data about 0.01-0.02 mag 
fainter according to theoretical models (Liu et al. 2000). 
Still, our measurements of the fluctuation apparent mag- 
nitudes tend to be slightly fainter, in part because of our 
different adopted parameters for the background galaxy 
LF (§ 3.2.2). 

To compute the /Tg-band SBF absolute magnitudes and 
SBF colors, we used I c data and /-band distances from 
Tonry et al. (2000b). The latter are calibrated using HST 
Ccpheid distances from Ferrarese et al. (2000a) for six 
nearby spiral galaxies which also have /-SBF measure- 
ments. The tabulated errors for the absolute /Tg-band 
SBF magnitudes include the errors from the /-SBF dis- 
tances added in quadrature. However, the listed errors do 
not include the systematic error in the HST Cepheid dis- 
tance scale, believed to be about ±0.16 mag (Mould et al. 
2000a) , which is mostly due to uncertainty in the distance 
to the Large Magellanic Cloud. 

We also computed Mj- S 's by using an HST Cepheid 
distance to the Fornax cluster. We adopted a Cepheid 
distance modulus of 31.42 ± 0.06 mag (19.2 Mpc) from 
the weighted average of distances to NGC 1326A and 
NGC 1365 (Ferrarese et al. 2000a); like Tonry et al. 
(2000a), we excluded NGC 1425 as part of Fornax. (We 
discuss the effect of the new Cepheid distances from Freed- 
man et al. 2001 in § 5.1.) We also added in quadrature 
a random error to account for the cluster's depth along 
the line of sight (0.06 mag). Just like the /-SBF cali- 
brated Mjf s 's, we did not include the systematic error 
in the HST Cepheid distance scale. Note that the total 
errors for the Cepheid-calibrated Mk s are smaller than 
the /-band SBF calibrated ones. This is because the er- 
ror added in quadrature from the Cepheid distance was 
only 0.08 mag, whereas the median error for the /-SBF 
distances to individual galaxies was 0.21 mag. 

If the Cepheid cluster distances are used, we are assum- 
ing the spiral galaxies used for the Cepheid measurements 
and the elliptical/lenticular galaxies with SBF data lie at 
the same distance. However, there is evidence that spirals 
in groups tend to lie at larger radii, with the ellipticals in 
the inner regions (Jacoby et al. 1992; Kelson et al. 2000). 
We circumvented this problem by choosing to use /-band 
SBF distances to individual galaxies. Although the /-SBF 
calibration is based on spirals with both HST Cepheid and 
/-SBF measurements, there are still some lingering con- 
cerns: SBF measurements in spiral bulges are challenging, 
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and bulges and ellipticals might have different stellar pop- 
ulations (e.g., Wyse et al. 1997). 

5. ANALYSIS 

Our new Kg-band sample of early-type galaxies in For- 
nax has more than doubled the number of high-quality IR 
SBF measurements. The sample covers an expanded range 
of galaxy properties (e.g., color, velocity dispersion, lumi- 
nosity, and morphological type) compared to past stud- 
ies. In particular, our new measurements include galaxies 
which are fainter and have bluer integrated colors. Also, 
all previous high-quality K-band SBF measurements were 
of ellipticals (except for the bulge of M31); we have added 
several SO galaxies. 

In the analysis which follows, we combine our Fornax 
sample with 10 galaxies in nearby clusters from the lit- 
erature. We use data for M32 and the bulge of M31 
from Luppino & Tonry (1993). Measurements for one Eri- 
danus cluster (NGC 1407) and five Virgo cluster galaxies 
(NGC 4365, 4406, 4472, 4552, 4636) come from Jensen 
et al. (1998). We do not include the lower S/N data of 
Jensen et al. given the concerns about systematic biases 
in SBF measurements at low S/N (see Jensen et al. 1996). 
For the same reason, we do not include data from Pahre & 
Mould (1994), and also because these authors did not cor- 
rect their SBF measurements for globular cluster contam- 
ination. Likewise, we hold our own Fornax data to similar 
standards and exclude mj- s measurements with S/N<3, 
namely those for IC 1919, NGC 1366, and NGC 1373. 

These literature data were all taken with the K' filter 
(1.9-2.3 /im; Wainscoat & Cowie 1992). To compare these 
to the slightly redder bandpass of our K$ filter, we sub- 
tracted 0.02 mag from the published data as suggested by 
the SBF models of Liu et al. (2000). We also changed 
the extinction corrections from the H I-dcrivcd values of 
Burstein & Heiles (1984) to the DIRBE/IRAS-derived val- 
ues of Schlegel et al. (1998) and used the latest V— I c 
galaxy colors from Tonry et al. (2000b). 

Finally, we include Ks-band data for NGC 3379 in Leo 
and NGC 4489 in Virgo from Mei et al. (2001). 

5.1. Calibration of Mk s 
5.1.1. Calculations 

Table 3 summarizes our various calibrations for M k s ■ 
There are several choices. We provide calibrations based 
on only our Fornax data or the total sample. Using only 
our Fornax data ensures the calibrations are based on mea- 
surements (and their associated errors) determined in a 
homogeneous fashion. For distances to galaxies, we use 
cither individual distances to galaxies from /-band SBF 
(Tonry et al. 2000b) or Ccphcid distances to the galaxy 
clusters as a whole. We adopt Ccpheid distance moduli 
of 24.44 ± 0.10 mag for M31 and M32, 30.08 ± 0.06 mag 
for Leo, and 31.03 ± 0.06 mag for Virgo (Ferrarese et al. 
2000b; see below for a discussion of the new distances by 
Freedman et al. 2001). 

In averaging the data, we either take a weighted average 
to derive a universal zero point for M x s , or we fit for a 

7 We test our algorithm using simulated data sets obeying y = A + Bx, where the errors in x and y are partially correlated. It is straight-forward 
to compute the shape of the tilted error ellipses (Cowan 1998). When the covariance is zero, our approach produces very similar results to 
the heuristic method of Press et al. (1992). For non-zero covariance, our approach is far more effective, producing essentially unbiased fits and 
accurate error estimates for A and B. Our simulations span a wide range of the parameter space in the coefficients, the sizes of the errors, and 
the amount of covariance, including values comparable to the observations. 



linear trend with V— I c to account for stellar population 
differences between galaxies, as is done for /-band SBFs. 
The data clearly favor a dependence on V — I c . A non- 
parametric evaluation using the Spearman rank-order cor- 
relation coefficient (r s ) finds a probability of only 3 x 10 -3 
for the dependence to occur by chance. Also, adopting a 
linear relation gives a lower reduced chi-square (x 2 ) than 
a universal zero point, and the derived slopes are statisti- 
cally different from zero. Finally, the existence of a color 
dependence is to be expected, given that /-band SBFs 
demonstrate substantial stellar population variations be- 
tween galaxies. 

We fit Mk s versus V— I c accounting for errors in both 
variables. For calibrations using the Ccpheid distances, we 
use the algorithm described in Press et al. (1992). How- 
ever, for calibrations using the /-band SBF distances, a 
different procedure is required. Since the /-band SBF 
distances depend on the V— I c color, or more specifically 
4.5 x (V— / c ), errors in V — I c will directly lead to errors 
in the resulting Mk s - We need to take into account this 
non-zero covariance when fitting for the relation between 
V — I c and Mk s - Such a covariance would bias the fitting 
towards a slope of 4.5, with more severe bias as the frac- 
tional errors in V — I c become larger. This effect was not 
accounted for by Jensen et al. (1998) in their fits. 

To avoid this bias, we start with the maximum likeli- 
hood approach described by Stetson (1992) which handles 
independent errors in both variables. We then enhance the 
method to account for the fact that the error ellipses for 
each datum are tilted, i.e., the covariance is non-zero. We 
validate our algorithm by applying it to simulated data 
sets. 7 These numerical experiments demonstrate that if 
we simply use an ordinary least-squares approach, or the 
method given by Press et al., the bias in the fits of M k s 
versus V—I c would small but not insignificant. The rea- 
son why the bias is relatively small is because of the high 
accuracy of the V — I c measurements. 

For all the calibrations, we used a standard cr-clipping 
scheme to avoid the undue influence of a few outliers. 
This resulted in the exclusion of NGC 1389 and 1419, 
both of which have very bright Ks-band fluctuations 
{Mk s ~ ~6.6 mag); these galaxies are discussed in the 
next section. More sophisticated iterative fitting schemes, 
e.g., based on robust weighting schemes described in Press 
et al., gave similar results. The main advantage of the cr- 
clipping scheme is its conceptual simplicity. 

Our preferred calibration relies on /-band SBF distances 
to individual galaxies, as described in § 4. Using the total 
sample with these distances gives 

M Ks = (-5.84±0.04) + (3.6±0.8)[(V-/ C ) -1.15] . (12) 

Figure 3 illustrates this calibration. Note that without ac- 
counting for the non-zero covariance between V — I c and 
Mk S i the fitted slope would have been 4.1 ± 0.8, with lit- 
tle change in the zeropoint. (As expected, if we do not 
account for the covariance, the fitted slope is biased to- 
wards a slope of 4.5.) There is a hint in the data that 



Liu, Graham, & Chariot 



the later-type galaxies may obey a different calibration. A 
fit to the 7 objects (6 SO's and the bulge of M3f ) finds a 
slope of 1.8±1.5. However, the sample is small and the 
significance is low. A larger sample is needed to test this 
possibility. 

5.1.2. Implications 

The dependence of M k s with V— I c color is a new find- 
ing. A hint of this effect was seen in the data set of 
Jensen et al. (1998), though it covered only about half 
the range in V — I c that our sample does. The depen- 
dence is seen clearly with the addition of Fornax galaxies 
with V — I c < 1.15. The Mk s slope is nearly as steep 
as for /-band SBFs (4.5 x (V-J c ); Tonry et al. 2000a), 
meaning accurate optical photometry improves the pre- 
cision of X-SBF distances. Moreover, the discovery of 
the slope eliminates published concerns about anomalously 
bright X-band SBFs for M32 (Luppino & Tonry 1993) 
and NGC 4489 in Virgo (Pahre & Mould 1994; Jensen 
et al. 1996; Mei et al. 2001). Both galaxies lie on the plot- 
ted calibration. (NGC 4489 has the bluest V—I c color 
in the sample, but if we remove it, the fitted zero point 
and slope of change negligibly.) However, we do find a 
few Fornax galaxies seem to have very bright Ks-h&nd 
fluctuations, more than 0.5 mag above the observed mean 
relation. We address the origin of these in the next section. 
Most of these galaxies are bluer and much less luminous 
than galaxies used for cosmological SBF distance measure- 
ments. 

The sense of the observed trend, with redder galaxies 
having fainter Ks-hemd fluctuations, is opposite that pre- 
dicted from theoretical models based on single-burst stel- 
lar populations. For the published SBF models, theoreti- 
cal calibrations based on ensemble averages of single-burst 
models (e.g., Worthey 1993) tend to produce brighter fluc- 
tuations for redder galaxies. The disagreement between 
these theoretical calibrations and the observations is not 
surprising given the relatively simply approach of averag- 
ing models with a range of ages and metallicites. Recently, 
Blakcslec et al. (2001) have generated more complicated 
models with three bursts of star formation. They select 
the ages and metallicities of the bursts to be consistent 
with the observed optical/IR SBF colors and the /-band 
SBF slope. Their resulting models have IR SBF magni- 
tudes which become brighter with redder integrated col- 
ors, with a slope of 2.9 in the /("-band, consistent with our 
observations. The zero point from their models is about 
~0.5 mag fainter than the observed calibration. (See § 5.2 
for more discussion about the theoretical models.) 

The discovery of a dependence of Mk s on V — I c color 
has important implications for Hq measurements using IR 
SBFs (Jensen et al. 1999; Jensen et al. 2001; Liu & Gra- 
ham 2001). In particular, Jensen et al. (2001) have used 
HST NICMOS F160W (1.6 fxm) data to measure SBFs in 
a sample of distant galaxies. They assumed the F160W 
SBFs were invariant from galaxy to galaxy. However, their 
set of nearby calibrators were on average ~0.05 mag bluer 
in V — I c than the set of distant galaxies used to deter- 
mine Hq- If the color dependence of F160W SBFs is as 
steep as that for Ks-b&nd SBFs, this means a difference 
of ^0.2 mag in the SBF absolute magnitudes between the 
two sets. The consequence would be that the measured 



distances for their distant set are overestimated by ~10% 
and likewise the resulting Hq underestimated. 

Calibrations of M k s based on Cepheid group distances 
have higher x 2 's since the Cepheid-calibrated M k s 's pos- 
sess much lower errors than /-SBF calibrated M^ s 's, as 
described in § 4. We have already added in quadrature 
an error term to account for the cluster depth (0.06 mag 
for Fornax and 0.08 mag for Virgo). To bring the Cepheid 
calibration of M k s to x 2 w 1, we would need to add in 
quadrature «0.15 mag of error. This result provides an 
estimate on the intrinsic cosmic scatter in M x s ■ (This 
calculation depends on the errors of rn~K s being accurate; 
these should be confirmed by re-observing at least a por- 
tion of the sample.) For comparison, the estimated cosmic 
scatter in /-band SBFs is about 0.05 mag (Tonry et al. 
1997), about a factor of three smaller. This difference is 
not surprising; stellar population models predict /-band 
SBFs will be degenerate in age and metallicity once cor- 
rected for differences in galaxy colors while /C-band SBFs 
will vary more strongly depending on the underlying pop- 
ulation (Worthey 1993; Liu et al. 2000). 

As mentioned before, this calibration does not account 
for any systematic error in the current HST Cepheid 
distance scale. For instance, recent comparisons of the 
Cepheid and maser distance measurements to the nearby 
spiral NGC 4258 suggest the Cepheid distances might be 
overestimated (Herrnstein et al. 1999; Maoz et al. 1999; 
Newman et al. 2001). In this work, we have used the 
Cepheid distances from Ferrarese et al. (2000a) since these 
are used to calibrate the /-band SBF distances compiled 
by Tonry et al. (2000b). Frccdman et al. (2001) have pre- 
sented new Cepheid distances by the HST Ho Key Project, 
which incorporate several changes in the analysis. These 
include use of an improved Cepheid period-luminosity re- 
lation, revised photometric calibrations, and a correction 
for the metallicity dependence of the Cephcids. The sign 
and amplitude of the metallicity correction remains quite 
uncertain. Using the new distances without the metal- 
licity correction leads to the /-SBF zero point becoming 
fainter by «0. 10-0. 15 mag. Our /tT,g-band SBF calibra- 
tion would change by the same amount. However with the 
metallicity correction adopted by Freedman et al., the I- 
SBF zero point is only slightly changed, becoming fainter 
by 0.06 mag. 

5.2. Recent Star Formation in Cluster Ellipticals 

The distribution of galaxies in Figure 3 has at least three 
characteristics: the elliptical and lenticular galaxies ap- 
pear to have comparable Mk s \ the JVs-band SBFs grow 
brighter as the integrated galaxy colors become bluer; and 
a few galaxies have SBFs much brighter than the bulk of 
the sample. These trends are also seen in Figure 4 which 
plots the Fornax cluster apparent SBF magnitudes, which 
are independent of the chosen galaxy distances, as a func- 
tion of galaxy properties. These phenomena provide new 
clues to the star formation histories of the cluster galaxies. 
While a complete analysis is deferred to a subsequent pa- 
per (M. Liu et al., in preparation), we address here some 
basic considerations. 

Since early-type galaxies follow a color-magnitude re- 
lation (more luminous ones have redder integrated colors 
[Visvanathan & Sandage 1977; Bower et al. 1992]), the cor- 
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relation between Mk s and V — I c color can also be seen 
as one between Mk s and galaxy luminosity/mass. Thus 
the bluer, less massive galaxies tend to have brighter K- 
band SBFs. As we discuss below, stellar population mod- 
els imply this trend reflects an age spread, with less mas- 
sive galaxies have more extended star formation histories. 
Furthermore, the bluer, lower luminosity galaxies (fainter 
than Mb ~ —19) span a wider range in Mk s - This sug- 
gests the star formation histories of these galaxies, in addi- 
tion to being more extended, are also more heterogeneous 
than those of the redder, more massive galaxies. A sim- 
ilar phenomenon may exist among early-type Virgo clus- 
ter galaxies, as the lower-mass [a < 100 km s^ 1 ) galax- 
ies exhibit much larger scatter in their Balmer and metal 
absorption lines than the more massive ones (Concannon 
et al. 2000). 

Figure 5 compares the Ks-h&nd SBFs with three sets 
of stellar population synthesis models, those of Liu et al. 
(2000), Worthey (1993), and Blakeslee et al. (2001). For 
the first set, the models used here are a slightly revised 
version of those published in Liu et al. -- see the Ap- 
pendix for details. The latter two models are calculated 
for slightly different if -band filters than used for our obser- 
vations, but the differences arc negligible for our purposes 
here. All the plotted models are single-burst stellar pop- 
ulations (SSPs), where the stars form coevally and then 
evolve passively. The models all use differing sets of evo- 
lutionary tracks and stellar spectral libraries (see discus- 
sions in Chariot et al. 1996, Liu et al. 2000, and Blakeslee 
et al. 2001). As a consequence, the SSP-equivalent ages 
and metallicitics inferred for the galaxies depend on the 
choice of models. 

The Liu et al. models indicate that the stellar popula- 
tions dominating the if -band SBFs of most galaxies are 
around solar metallicity with a factor of 3-4 spread in 
age, including some as young as «2 Gyr. Given the mea- 
surement errors in M k s , the data also allow for about a 
factor of two spread in metallicity. The Worthey mod- 
els offer qualitatively similar results, though with a lower 
mean metallicity. Unfortunately, the Worthey models are 
unavailable for young, metal-poor populations so they do 
not cover the full observational locus. The Blakeslee et al. 
models suggest a comparable spread in age and a much 
wider spread in metallicity, extending above their maxi- 
mum available metallicity (1.5x solar). These models sug- 
gest the effects of age and metallicity are nearly degenerate 
for metal-rich SSPs. It is worth noting that despite the dis- 
agreements in the model predictions, all three sets of mod- 
els in Figure 5 indicate ages younger than 5 Gyr for many 
galaxies. Such young ages imply these galaxies have very 
luminous giant branches, as is seen in intermediate-age 
(«0.5-5 Gyr) Magellanic Cloud clusters (Mould & Aaron- 
son 1979; Frogel et al. 1980; Mould & Aaronson 1980; Co- 
hen et al. 1981). 

The disagreement is most severe between the Liu et al. 
and Blakeslee et al. predictions, which are derived from 
the Bruzual & Chariot (2000) and Vazdekis et al. (1996) 
models, respectively. Note that the two sets of models 
are computed for different metallicities; in particular, the 
Blakeslee et al. models are not as metal-rich. At metal- 
licities of solar and above, which are most relevant to el- 
liptical galaxies, the Liu et al. models predict the if-band 



SBFs grow brighter with increasing metallicity, as would 
be expected based on the steady increase with metallic- 
ity in the if-band magnitude of the tip of the red giant 
branch for Galactic globular clusters (Ferraro et al. 2000). 
In contrast, Blakeslee et al. predict the IR SBFs do not 
change much between solar and super-solar metallicity. 
For the reddest ellipticals (V — I c > 1.2), the metallici- 
ties inferred in Figure 5 from their models are very large. 
Furthermore, for SSPs younger than «8 Gyr, Blakeslee et 
al. predict that the IR SBF magnitudes decrease with in- 
creasing metallicity. This "inversion" is surprising given 
two known aspects of AGB stars, which are expected to 
contribute a substantial fraction of the IR SBF signal (see 
Liu et al. 2000). First, the main-sequence turnoff mass in- 
creases with metallicity at fixed age (Renzini & Fusi Pecci 
1988). Since AGB stars follow a core mass- luminosity rela- 
tion, higher metallicity is expected to lead to brighter AGB 
stars. Likewise, solar-metallicity intermediate-age clusters 
are predicted to have AGB stars with higher luminosities 
than lower metallicity clusters (see Silva & Bothun 1998 
and Liu et al. 2000). Both of these argue for brighter AGB 
stars with increasing metallicity, hence brighter infrared 
M's. For these reasons, we favor the Liu et al. models, 
though the issue remains to be resolved. _ 

Figure 6 compares the SSP models with I c —Ks SBF 
color measurements, which are independent of the galaxy 
distances. Qualitatively the results are the same: the 
Liu et al. and Worthey models give a large spread in age 
with roughly solar metallicity (as in the case for Mk s -> 
the errors in the SBF colors allow for some spread in 
the metallicity); the Blakeslee et al. models find a com- 
parable spread in age with a larger spread in metallicity. 
These comparisons using l c -Ks lead to slightly larger age 
spreads and slightly smaller metallicity spreads than if we 
use M k s ■ The fluctuation colors for all the models follow 
the same trend, becoming redder primarily as the metal- 
licity increases. The Blakeslee et al. models for I c -Ks are 
more similar to the other models and better behaved than 
for Mk s \ this further suggests there may be a problem 
with their IR fluctuation magnitudes. 

Kuntschner (2000) has recently studied the star forma- 
tion history of early-type galaxies in the Fornax cluster 
using optical absorption lines. He uses Balmer and metal 
lines to infer galaxian ages and metallicities, respectively. 
(See Trager et al. 2000 for an independent analysis of the 
same data set.) There are 13 galaxies in common between 
our work and his, most of which are ellipticals and not 
SO's. The faintest SO's in his sample, which show a very 
large spread in age, are not in our SBF sample. The size 
of the age and metallicity spreads inferred from his ab- 
sorption line analysis and our SBF data are in reasonable 
agreement. He infers a range in [Fe/H] of -0.25 to +0.25, 
compatible with our data especially given its errors. He 
finds an age range of 5-12 Gyr, which is a slightly smaller 
spread and a higher mean age than our results neither 
of these differences is surprising given that the SBF data 
are more sensitive to recent star formation (see below). 

The above comparisons use single-burst population 
models. But since the star formation history of early-type 
galaxies remains uncertain, it is important to consider the 
effect of multiple episodes of star formation. In low red- 
shift clusters, these galaxies follow a color-magnitude re- 
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lation with little intrinsic scatter, which is typically in- 
terpreted as a small spread in age (Bower et al. 1992). 
However, intermediate-age (few Gyr) stellar populations 
can exist in small amounts provided the bulk of the stars 
are old (Bower et al. 1998). Indeed, there are several lines 
of evidence from distant (z > 0.3) clusters that early-type 
galaxies have experienced more than a single episode of 
star formation (e.g., Chariot & Silk 1994; Barger et al. 
1996; Poggianti et al. 1999; Ferreras & Silk 2000), with 
the majority of the stellar mass forming at high redshift 
but a minority fraction originating more recently. 

This evidence is based on optical colors and absorp- 
tion line measurements, which have several intrinsic dif- 
ferences from SBF data for stellar populations. Absorp- 
tion lines are typically measured only in the central re- 
gions of galaxies, while SBF measurements sample a much 
larger volume. Also, the light from main-sequence stars 
contribute to the optical colors and absorption lines (e.g., 
Bruzual & Chariot 1993), but the optical/IR SBF signal is 
much more dependent on the giant stars. Because of this 
weighting to luminous cool stars, changes in metallicity 
might be more easily discerned from SBF measurements 
than from integrated spectral properties (Liu et al. 2000). 
Also, the difference in weighting means that model predic- 
tions for optical colors/lines and SBFs are subject to very 
different uncertainties. For instance, the interpretation of 
metal absorption lines is hampered by the enhancement 
of a-element in ellipticals relative to the models (Worthey 
et al. 1992) and the limited range in stellar temperatures 
and gravities used for model inputs (e.g., Worthey et al. 
1994; Maraston & Thomas 2000), neither of which have 
a direct bearing on SBF predictions. Hence, concordance 
between these lines of evidence would help confirm the 
robustness of the results. 

In a scenario where most of the stellar mass forms at 
high redshift with a small fraction forming more recently, 
we would expect the following. For recent star formation 
occurring with the same metallicity as the older in situ 
population, galaxies will move in the {V— I c , Mk s } plane 
along lines of constant metallicity. Thus, the spread in 
age inferred from the Liu et al. and Worthey models could 
arise from late bursts of star formation with a common 
metallicity. On the other hand, if a late burst occurs with 
a different metallicity, we expect noticeable changes to the 
position of the galaxies in the {V— I c , Mk s } plane. De- 
spite its relatively small mass, the younger population will 
be a substantial contributor to the total luminosity for 
the first few Gyr because of its much higher light-to-mass 
ratio. Since SBFs are sensitive to the second moment of 
the stellar luminosity function, we expect the SBFs to be 
dominated by this younger population. Ages inferred from 
SBFs would be expected to be younger and have larger 
scatter than those from integrated colors/lines. Models 
predict that the IR SBFs depend on metallicity (Figures 5 
and 6), albeit to varying degrees. The Worthey and Liu et 
al. models predict a strong metallicity dependence at all 
ages; if correct, then in the first few Gyr after the burst, 
the IR SBFs could be either brighter or fainter, depending 
on the metallicity of the burst relative to the underlying 
population. 

Figure 7 illustrates the effect of recent star formation 
using evolving population models of Bruzual & Chariot 



(2000). The majority of the model population has solar 
metallicity. A second burst of star formation is assumed 
to occur 6 Gyr after formation and involving 20% of the 
final galaxy mass. When the two bursts have the same 
metallicity, the effect of the second burst is to reduce the 
model-inferred age with little effect on the inferred metal- 
licity. However, when the second burst has a higher metal- 
licity (in this case [Fe/H]=+0.4), M k s is greatly bright- 
ened while the range in V — I c galaxy color is unchanged. 
Notice that the effect of the second burst on Mk s lingers 
for many Gyr. This is in sharp contrast to ordinary op- 
tical colors, which revert to their pre-burst state in only 
a few Gyr (e.g., Chariot & Silk 1994). Similarly, for a 
second burst with a much lower metallicity ([Fe/H]=— 1.7, 
typical of local dwarf galaxies), the range in optical col- 
ors is similar to the other models only ~2 Gyr after the 
burst, but the IR SBFs become much fainter for long af- 
terwards. Therefore, the offset in SBF magnitudes/colors 
might be able to distinguish late bursts of star formation 
which have very different metallicities from the main pop- 
ulation, even after the perturbations to the optical colors 
and absorption lines have subsided. 

The two-burst models shown in Figure 7 should be taken 
as illustrative, not definitive. They do show that if a small 
fraction of the total galaxy mass was formed in the last few 
Gyr, this would suffice to explain the spread of the opti- 
cal/IR SBFs and integrated colors. In addition, the mod- 
eling suggests that the galaxies with unusually bright K$- 
band SBFs, such as NGC 1389 and 1419, result from recent 
star formation involving gas which was enriched in metals. 
Likewise, the absence of galaxies with unusually faint Kg- 
band SBFs implies a lack of metal-poor intermediate-age 
stars, which would form from unenriched infalling gas. A 
more complete analysis of the stellar populations of these 
galaxies is the subject of a future paper. 

6. CONCLUSIONS 

We have presented ifg-band SBF data for 19 early- type 
Fornax cluster galaxies. This doubles the number of high- 
quality IR SBF measurements to date. Combining our 
measurements with data from the literature for ellipticals 
in nearby clusters, we have calibrated M k s as a distance 
indicator. We offer calibrations using either HST Cepheid 
distances to the clusters or /-band SBF distances to indi- 
vidual galaxies. When using the latter, we account for the 
covariance between V — I c and Mk s - For both options, 
any systematic change in the HST Cepheid distance scale 
will directly impact the M k s calibration. 

We have found /Cg-band fluctuation magnitudes vary 
considerably between galaxies. However, we also find that 
the Ks-band SBFs correlate with V — I c galaxy color, like 
/-band SBFs. The existence of this correlation means Kg- 
band distances can be improved by optical color measure- 
ments to correct for stellar population variations between 
galaxies. It also suggests that the Hq measurement by 
Jensen et al. (2001) using IR SBFs might be underesti- 
mated by ~10% since they assumed no color-dependence 
for the SBF magnitudes. The later-type galaxies may fol- 
low a different correlation than the ellipticals, though this 
is uncertain with the existing data. In addition, this find- 
ing resolves published concerns that the SBFs of NGC 4489 
and M32 have anomalously bright /f-band SBFs: their 
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brighter fluctuations are in accord with their bluer inte- 
grated colors. Overall, the intrinsic scatter in /T-band SBF 
appears to be significantly larger than /-band SBF. 

We also find a few galaxies have very bright /Cs-band 
fluctuations, more than 0.5 mag above the observed mean 
relation. Most of these galaxies are bluer and much less lu- 
minous than the bright early-type galaxies used for cosmo- 
logical distance measurements. Nevertheless the existence 
of such galaxies may be a concern for attempts to measure 
distances and determine Hq with IR SBFs. Larger sam- 
ples are needed to determine the prevalence and extent of 
this phenomenon. 

The existence of a correlation in M k s (and I c -Ks) with 
V — I c galaxy color is a new clue into the star formation 
histories of cluster galaxies. Such a relation is known to ex- 
ist for /-band SBFs. However, stellar population models 
predict the effects of age and metallicity are largely de- 
generate in this bandpass, and hence /-band SBFs alone 
are expected to be of little use for stellar population stud- 
ies. On the other hand, models predict strong effects from 
stellar population variations on IR fluctuation magnitudes 
and optical/IR fluctuation colors. Hence, our discovery of 
a systematic relation between IR SBFs and galaxy color 
means we can determine specific age/metallicity combina- 
tions, which then will be a reflection on the star formation 
histories. Given the current uncertainties in modeling the 
RGB and AGB stars which dominate the IR SBF signal, 
these interpretations will depend on the choice of stellar 
population models. We argue the SBF model predictions 
from Liu et al. (2000) are probably most reasonable based 
on what is known from resolved stellar populations in Lo- 
cal Group star clusters. We also point out some potential 
shortcomings with the models of Blakeslee et al. (2001). 

Both the Liu et al. and Worthey (1994) models suggest 
that the trend in Mk s with V — I c originates from varia- 
tions in the age of the populations dominating the SBF sig- 
nal. Most of the galaxies are inferred to have roughly solar 
metallicity, with perhaps a spread of a factor of two. Simi- 
lar conclusions are reached using I c ~Ks fluctuation colors, 
which are independent of galaxy distances. All the mod- 
els suggest the single-burst equivalent ages for the bluest 
galaxies are intermediate-age (<5 Gyr), implying these 
galaxies have very luminous extended giant branches, sim- 
ilar to those found in intermediate-age Magellanic Cloud 
clusters. In the context of scenarios where star formation 
occurs in more than one episode, the SBF/galaxy color 
trend can be explained by the occurrence of late bursts 
of star formation with metallicity similar to the older in 
situ population comprising most of the stellar mass. Such 
bursts act to change the optical colors, which reduces the 
inferred age, without causing the fluctuations to deviate 
from the main trend. 

In a similar fashion, the few galaxies with very bright IR 
SBFs might be explained by late bursts of star formation 
with a higher metallicity than the old stars. In this pic- 



ture, the lack of galaxies with much fainter IR fluctuations 
disfavors secondary star formation from more metal-poor 
gas. The existence of galaxies with very bright fluctu- 
ations among the bluer, less luminous galaxies indicates 
that the star formation histories of these galaxies were 
more heterogeneous than those of the redder, more massive 
galaxies. This finding may seem at first to be at variance 
with scenarios of hierarchical galaxy formation in a cold 
dark matter dominated cosmology (e.g. Baugh et al. 1996; 
Kauffmann & Chariot 1998), which predict that the more 
massive galaxies have more extended formation histories. 
However, in hierarchical cosmologies merging rates depend 
strongly on mass — low-mass galaxies experience fewer 
recent mergers than high-mass galaxies. If mergers trig- 
ger efficient star formation, Kauffmann et al. (2001) find 
that lower-mass galaxies will experience much more widely 
spaced bursts than massive galaxies, and hence their star 
formation histories will be more heterogenous. 

We caution that our results for the calibration of K$- 
band SBFs and the inferred stellar populations are for 
galaxies in nearby clusters. It remains an open question 
how these results depend on environment. Since IR SBFs 
are predicted by all models to strongly depend on age 
and metallicity, they could vary substantially with envi- 
ronment, though the data for Local Group galaxies seem 
to agree well with those for galaxies in the much denser 
Virgo and Fornax clusters. Expanded samples will be valu- 
able both for strengthening the calibration of IR SBFs as 
cosmological distance indicators and for deciphering the 
stellar population histories of early-type galaxies. 
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APPENDIX 

REVISED SBF MODELS 



In the version of the Bruzual-Charlot models used here, the thermally pulsing phase of the asymptotic giant branch 
(TP-AGB) has been slightly improved relative to the version used in Liu et al. (2000). The models rely on the same 
calculations and spectral calibrations of TP-AGB stars (M-type, C-type, and superwind phase) as described by Liu et al. 
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However, the sampling of these stars has been improved to better account for the extended giant branches of Magellanic 
Cloud and Galactic Bulge clusters, as observed by Ferraro et al. (1995) and Guarnieri et al. (1998). Noticeable changes 
occur in all the SBF magnitudes for the young (1-3 Gyr) models with sub-solar metallicities and also in the optical SBFs 
(V R c I c ) for the metal-rich (Z > Z Q ) models. Table Al presents the new set of default models, which supersedes Table 2 
of Liu et al. 
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Fig. 1. — Top: Ks-haxid images of two Fornax galaxies, each 1.5' on a side. Bottom: Same images after subtracting an 
elliptical model. The innermost region has been masked. SBFs are seen as the faint mottling near the center, with a 
larger amplitude for higher surface brightness regions. Globular clusters and background galaxies are also present. 
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Fig. 2. — Ks-band fluctuation power spectra for the central 256x256 pixel region of our Fornax galaxy images. The 
galaxy power spectra are fitted by the sum (solid line) of a scaled version of the PSF (P x E(k)) and a constant 
(Pi); the dashed lines show the contributions from these two components. We use wavenumbers from k = 20 (spatial 
frequencies of 0.05 pixel -1 ) to k = 128 (the Nyquist frequency) and do the fit with the two-dimensional power spectra. 
The one-dimensional azimuthal averages are plotted. 
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Fig. 3. — Combined M k s data set for galaxies with high-quality SBF measurements plotted against integrated V— I c 
galaxy colors. Distances to the individual galaxies come from /-band SBF (Tonry et al. 2000b). The open symbols are 
SO's and spiral bulges (M31), and the filled ones are ellipticals. The Fornax data are from this work; the Local Group 
(M31 and M32, with M32 being the galaxy with bluer V—I c ) data are from Luppino & Tonry (1993); and the Virgo, 
Eridanus, and Leo data are from Jensen et al. (1998) and Mei et al. (2001). The values in parentheses are the number of 
galaxies in each group with SBF measurements. The dashed line is the best fit. The fit excludes the two Fornax galaxies 
with very bright fluctuations (M Ks ~ -6.6 mag), NGC 1419 (E) and NGC 1389 (SO). The reddest galaxy (lower right 
of plot) is NGC 4636 in Virgo. 
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Fig. 4. — Ks-band SBF apparent magnitudes for our Fornax cluster sample as a function of galaxy properties. Top: 
V — I c galaxy color. Bottom: Absolute i?-band magnitude, with apparent magnitudes Bt from de Vaucouleurs et al. 
(1991) and distances and errors from /-band SBF data of Tonry et al. (2000b). Elliptical galaxies are listed in roman 
type and SO's in italics. The errors in mjf s for the brighter half of the sample (M(Bt) < —19) are comparable to the 
height of the labels; for the fainter half, the errors are a few to several times larger than the label height. Errors in 
M(Bt) are comparable to or smaller than the width of the labels. 
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Fig. 5. — Comparison of Ks-band SBF measurements with three sets of stellar population synthesis models. The data 
are the same as in Figure 3, using /-band SBF distances to individual galaxies. Filled circles (•) are ellipticals, and open 
circles (o) are lenticular galaxies and spiral bulges. The average error bars are plotted in the upper right corners. Models 
of a fixed metallicity have the same symbol, with increasing symbol size representing increasing age; solar metallicity 
models are plotted with a diamond (Z = 0.02, 0.02, and 0.019). Line connect models with the same age. The Liu et al. 
(2000) models are 2, 3, 5, 8, 12, and 17 Gyr old. The Worthey (1994) models cover the same, except they are unavailable 
for young, metal-poor populations. The models from Blakeslee et al. (2001) are 4, 5, 8, 13, and 18 Gyr; note that their 
metal-rich models are plotted with the same symbols as the other models but have different Z values. 
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Fig. 6. — Comparison of I c —Kg SBF colors with stellar population synthesis models. The data and models are the 
same as in Figure 5. Solar metallicity models are plotted with a diamond (Z = 0.02, 0.02, and 0.019). Filled circles (•) 
are ellipticals, and open circles (o) are lenticular galaxies and spiral bulges. The average error bars are plotted in the 
upper left corners. 
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Fig. 7. — An illustration of the effect of recent star formation on the i^s-band SBF magnitudes and integrated galaxy 
colors using the stellar population models of Liu et al. (2000). Filled circles (•) are data for ellipticals, and open circles (o) 
are for lenticular galaxies. The average error bars are plotted in the upper left corner. The solid line shows a single-burst 
solar-metallicity galaxy model aging from 2 to 12 Gyr, moving from left to right on the plot. The broken lines show 
models where a second burst of star formation occurs 6 Gyr after the first burst and amounts to 20% of the final mass. 
When the second burst has the same metallicity as the first burst (dot-dashed line), the effect of the burst is to shift 
the model- inferred age roughly along the same locus. When the second burst has a higher metallicity ([Fe/H]=+0.4; 
upper dashed line), the Mks lii S rea tly increased while the range of V — I c galaxy color is unaffected. When the second 
burst has a lower metallicity ([Fc/H]=— 1.7; lower dotted line), the Mk s is greatly reduced. The open squares are spaced 
1 Gyr apart, starting at 7 Gyr age (i.e., 1 Gyr after the burst occurs). The oldest (12 Gyr) version of each model is 
marked with a filled square. 
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For 
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Name 


Type a 


'T'b 


Bj* 


Dia. a 


a c 


A B d 


V-I c c 








(mag) 


(arcmin) 


(km s -1 ) 


(mag) 


(mag) 


IC 1919 


SA(rs)0- 


-3 


13.80 


1.6 


X 


1.2 


26 


0.06 


1.108 ± 0.018 


IC 2006 


E 


-4 


12.21 


2.1 


X 


1.8 


122 


0.05 


1.183 ± 0.018 


NGC 1336 


SA0- 


3 


13.10 


2.1 


X 


1.5 


115 


0.05 


1.124 ± 0.032 


NGC 1339 


E4 


-4 


12.51 


1.9 


X 


1.4 


171 


0.06 


1.134 ± 0.012 


NGC 1344 f 


E5 


-5 


11.27 


6.0 


X 


3.5 


159 


0.08 


1.135 ± 0.011 


NGC 1351 


SAO- pec 


3 


12.46 


2.8 


X 


1.7 


147 


0.06 


1.148 ± 0.016 


NGC 1366 


SO 


2 


12.86 


2.1 


X 


0.9 




0.07 


1.095 ± 0.018 


NGC 1373 


E 


-4 


14.12 


1.1 


X 


0.9 


79 


0.06 


1.085 ± 0.013 


NGC 1374 


E 


-4 


12.00 


2.5 


X 


2.3 


207 


0.06 


1.146 ± 0.016 


NGC 1375 


SABO 


2 


13.18 


2.2 


X 


0.9 


69 


0.06 


1.070 ± 0.019 


NGC 1379 


E 


5 


11.80 


2.4 


X 


2.3 


130 


0.05 


1.143 ± 0.019 


NGC 1380 


SAO 


2 


10.87 


4.8 


X 


2.3 


240 


0.08 


1.197 ± 0.019 


NGC 1380BS 


SAB(s)0- 


3 


13.87 


1.5 


X 


1.3 


99 


0.07 


1.106 ± 0.013 


NGC 1387 


SAB(s)0- 


3 


11.68 


2.8 


X 


2.8 




0.06 


1.208 ± 0.047 


NGC 1389 


SAB(s)0 


3 


12.42 


2.3 


X 


1.4 




0.05 


1.145 ± 0.019 


NGC 1399 


El pec 


-5 


10.55 


6.9 


X 


6.5 


359 


0.06 


1.227 ± 0.016 


NGC 1404 


El 


-5 


10.97 


3.3 


X 


3.0 


242 


0.05 


1.224 ± 0.016 


NGC 1419 


E 


-5 


13.48 


1.1 


X 


1.1 


129 


0.06 


1.110 ± 0.018 


NGC 1427 


E5 


-4 


11.77 


3.6 


X 


2.5 


170 


0.05 


1.152 ± 0.018 



a From the NASA/IPAC Extragalactic Database (NED). 

b From dc Vaucoulcurs ct al. (1991) (RC3). 

c Data are from Hypercat (http: //www-obs . vmiv-lyonl . f r/hypercat), which uses an up- 
dated version of the literature compilation by Prugniel & Simien (1996). The exception is 
NGC 1373; its datum is from Graham et al. (1998). 

d From Schlegel et al. (1998). 
°From Tonry et al. (2000b). 

'a.k.a. NGC 1340. 
Sa.k.a. NGC 1382. 
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Table 2 












Ks-BAND 


SBF Measurements 






Galaxy 


(tJ-Ks) 


(P -P r )/Pl 


m Ks 


M Ks 


-=— ?Cepheid 

M Ks P 


Ic-K s 




(mag arcsec~ 2 ) 




(mag) 


(mag) 


(mag) 


(mag) 


IC 1919 


18.52 


6.0 ± 2.5 


24.84 ± 0.50 


-6.47 ± 0.53 


-6.58 ± 0.50 


4.54 ± 0.52 


IC 2006 


17.77 


4.7 ± 0.3 


25.76 ± 0.11 


-5.83 ± 0.31 


-5.66 ± 0.14 


4.23 ± 0.29 


NGC 1336 


18.49 


2.9 ± 0.4 


25.71 ± 0.23 


-5.67 ± 0.31 


-5.71 ± 0.25 


3.82 ± 0.28 


NGC 1339 


17.57 


4.3 ± 0.4 


25.85 ± 0.08 


-5.76 ± 0.36 


-5.57 ± 0.12 


3.94 ± 0.36 


NGC 1344 


17.21 


18.6 ± 0.9 


25.59 ± 0.10 


-5.89 ± 0.31 


-5.83 ± 0.12 


4.08 ± 0.31 


NGC 1351 


17.60 


5.4 ± 0.6 


25.93 ± 0.11 


-5.68 ± 0.19 


-5.48 ± 0.13 


3.93 ± 0.17 


NGC 1366 


19.00 


2.2 ± 0.6 


25.51 ± 0.43 


-6.11 ± 0.52 


-5.91 ± 0.43 


4.12 ± 0.51 


NGC 1373 


18.56 


9.0 ± 1.9 


25.11 ± 0.37 


-6.67 ± 0.60 


-6.31 ± 0.38 


4.63 ± 0.59 


NGC 1374 


17.05 


8.3 ± 0.3 


25.66 ± 0.05 


-5.82 ± 0.14 


-5.75 ± 0.09 


4.06 ± 0.11 


NGC 1375 


17.82 


10.6 ± 1.2 


25.52 ± 0.26 


-6.06 ± 0.29 


-5.90 ± 0.27 


3.96 ± 0.28 


NGC 1379 


17.08 


12.0 ± 0.9 


25.66 ± 0.08 


-5.85 ± 0.17 


-5.75 ± 0.12 


4.08 ± 0.14 


NGC 1380 


16.28 


7.4 ± 0.5 


25.39 ± 0.04 


-5.84 ± 0.18 


-6.02 ± 0.09 


4.31 ± 0.15 


NGC 1380B 


17.95 


5.3 ± 0.4 


25.78 ± 0.25 


-5.99 ± 0.40 


-5.63 ± 0.27 


4.05 ± 0.39 


NGC 1387 


17.28 


8.2 ± 0.7 


25.65 ± 0.06 


-5.89 ± 0.27 


-5.77 ± 0.10 


4.41 ± 0.15 


NGC 1389 


17.53 


14.0 ± 0.8 


25.07 ± 0.10 


-6.61 ± 0.21 


-6.35 ± 0.13 


4.85 ± 0.19 


NGC 1399 


16.01 


5.2 ± 0.5 


26.14 ± 0.12 


-5.36 ± 0.20 


-5.28 ± 0.15 


3.97 ± 0.18 


NGC 1404 


16.67 


8.9 ± 0.6 


25.89 ± 0.05 


-5.72 ± 0.20 


-5.53 ± 0.10 


4.31 ± 0.17 


NGC 1419 


17.92 


16.6 ± 1.4 


24.65 ± 0.13 


-6.77 ± 0.27 


-6.77 ± 0.15 


4.85 ± 0.26 


NGC 1427 


17.51 


9.5 ± 0.5 


25.46 ± 0.06 


-6.40 ± 0.25 


-5.96 ± 0.10 


4.67 ± 0.23 



Note. -- The quantity (Po — Pr)/Pi measures the S/N of the stellar SBF signal (i.e., includes correction for 
unresolved contaminating point sources). The tabulated mx s has been corrected for extinction. Ks-b&nd M'b 
were computed using two distance measurements to the galaxies: J-band SBF distances to individual galaxies and 
a Cepheid distance to the cluster as a whole. The errors on the Mk s from the Cepheid group distance includes 
an allowance for the cluster's rms depth but does not include the systematic error in the HST Cepheid zero point, 
estimated to be ±0.16 mag (Mould et al. 2000a). 



Table 3 
Calibration of Mkc 



Clusters 



Distances 



N 



e 



Fornax only /-band SBF (individual) 
Cepheid (cluster) 

All galaxies /-band SBF (individual) 
Cepheid (cluster) 21 



-5.81 ± 0.06 




14 


0.23 


1.0 


-5.84 ± 0.06 


2.0 ± 1.4 


14 


0.21 


0.9 


-5.74 ± 0.03 




14 


0.20 


2.8 


-5.78 ± 0.04 


4.1 ± 1.6 


14 


0.22 


2.5 


-5.76 ± 0.04 




24 


0.28 


2.1 


-5.84 ± 0.04 


3.6 ± 0.8 


24 


0.22 


1.2 


-5.68 ± 0.03 




23 


0.25 


4.2 


-5.75 ± 0.03 


4.3 ± 0.7 


23 


0.23 


2.4 



^oes not include NGC 1407 in Eridanus since no Cepheid distance is available. 



Note. — Calibration of M k s using only our Fornax results and our results combined 
with revised published data for Virgo, Leo, Eridanus, and Local Group galaxies. We 
cither use /-band SBF distances to individual galaxies or mean distances to the clus- 
ters from Cepheid distances in Ferrarese et al. (2000a). The linear fits arc of the form 
M k s = a + b [(V — I c )o — 1.15], and the number of galaxies used for the fit is tabulated 
(N) . Also listed are the resulting rms (in mags) of the points after the fit and the reduced 
chi-square (x 2 ) ■ See § 5.1 for details. 
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Revised BC2000 Model Predictions 


























Gyr 


B 


V 


Re 


U 








j 




H 


K' 


K s 


K 




L 


U 


M 


V-I c 


V-K 


J-K 


H/3 


Mg 2 


Mgb 


Hja 


C4668 




z 


F814W 


F110M 


F110W 


F160W 


F222M 




0.0001 


1 


0.06 


-0.56 


-1.27 


-2.24 


-2.15 


-3.24 


-3.22 


-3.66 


-4.90 


-5.07 


-5.58 


-5.60 


-5.65 


-5.57 


-7.26 


-7.29 


-7.58 


0.41 


1.15 


0.43 


6.72 


-0.04 


-2.28 


8.79 


0.97 




0.0001 


2 


0.35 


-0.78 


-1.66 


-2.63 


-2.54 


-3.58 


-3.56 


-3.96 


-5.11 


-5.26 


-5.77 


-5.78 


-5.82 


-5.75 


-6.83 


-6.89 


-7.10 


0.69 


1.77 


0.60 


4.72 


-0.04 


-1.31 


6.51 


0.38 




0.0001 


3 


0.46 


-0.68 


-1.54 


-2.47 


-2.39 


-3.39 


-3.37 


-3.77 


-4.91 


-5.06 


-5.58 


-5.59 


-5.63 


-5.56 


-6.55 


-6.61 


-6.81 


0.73 


1.82 


0.60 


4.22 


-0.03 


-0.98 


6.07 


0.20 




0.0001 


5 


0.46 


-0.70 


-1.50 


-2.28 


-2.22 


-3.03 


-3.00 


-3.36 


-4.37 


-4.51 


-5.06 


-5.07 


-5.12 


-5.01 


-6.03 


-6.06 


-6.32 


0.78 


1.85 


0.57 


3.48 


-0.01 


-0.47 


4.86 


-0.08 




0.0001 


8 


0.46 


-0.66 


-1.44 


-2.14 


-2.09 


-2.69 


-2.66 


-2.93 


-3.58 


-3.66 


-3.93 


-3.94 


-4.00 


-3.95 


-5.28 


-5.29 


-5.60 


0.81 


1.86 


0.54 


3.10 


0.00 


-0.13 


3.92 


-0.24 




0.0001 


12 


0.50 


-0.58 


-1.33 


-2.00 


-1.95 


-2.52 


-2.49 


-2.76 


-3.38 


-3.46 


-3.70 


-3.71 


-3.75 


-3.72 


-4.89 


-4.89 


-5.18 


0.85 


1.93 


0.55 


2.91 


0.02 


0.09 


3.08 


-0.33 




0.0001 


17 


0.51 


-0.56 


-1.29 


-1.94 


-1.89 


-2.45 


-2.42 


-2.68 


-3.28 


-3.36 


-3.55 


-3.57 


-3.59 


-3.58 


-4.44 


-4.45 


-4.68 


0.89 


2.00 


0.57 


2.73 


0.03 


0.27 


2.38 


-0.42 




0.0004 


1 


0.41 


-0.63 


-1.45 


-2.33 


-2.26 


-3.21 


-3.18 


-3.57 


-4.65 


-4.80 


-5.28 


-5.29 


-5.34 


-5.27 


-7.00 


-7.04 


-7.34 


0.58 


1.47 


0.49 


5.50 


-0.02 


-1.24 


6.51 


-0.39 




0.0004 


2 


0.27 


-1.03 


-1.92 


-2.83 


-2.75 


-3.71 


-3.69 


-4.07 


-5.14 


-5.29 


-5.77 


-5.78 


-5.81 


-5.75 


-6.79 


-6.84 


-7.06 


0.77 


1.91 


0.62 


4.08 


-0.01 


-0.52 


5.09 


-0.47 




0.0004 


3 


0.47 


-0.82 


-1.68 


-2.58 


-2.50 


-3.45 


-3.43 


-3.82 


-4.90 


-5.04 


-5.54 


-5.55 


-5.59 


-5.52 


-6.51 


-6.56 


-6.78 


0.79 


1.93 


0.62 


3.75 


-0.00 


-0.30 


4.76 


-0.54 




0.0004 


5 


0.51 


-0.83 


-1.70 


-2.51 


-2.45 


-3.27 


-3.25 


-3.60 


-4.57 


-4.71 


-5.19 


-5.20 


-5.24 


-5.17 


-6.06 


-6.10 


-6.33 


0.84 


1.98 


0.61 


3.13 


0.01 


0.15 


3.27 


-0.58 




0.0004 


8 


0.53 


-0.79 


-1.60 


-2.30 


-2.25 


-2.88 


-2.85 


-3.13 


-3.80 


-3.89 


-4.14 


-4.15 


-4.19 


-4.16 


-5.34 


-5.35 


-5.64 


0.88 


2.00 


0.58 


2.66 


0.03 


0.50 


1.81 


-0.61 




0.0004 


12 


0.55 


-0.74 


-1.54 


-2.23 


-2.19 


-2.81 


-2.77 


-3.06 


-3.71 


-3.79 


-4.01 


-4.01 


-4.03 


-4.03 


-4.90 


-4.91 


-5.12 


0.91 


2.05 


0.59 


2.56 


0.04 


0.64 


1.53 


-0.64 




0.0004 


17 


0.55 


-0.72 


-1.52 


-2.21 


-2.16 


-2.78 


-2.75 


-3.03 


-3.68 


-3.77 


-3.96 


-3.96 


-3.97 


-3.98 


-4.65 


-4.66 


-4.80 


0.93 


2.10 


0.60 


2.64 


0.05 


0.66 


1.77 


-0.66 




0.0040 


1 


1.20 


-0.12 


-1.33 


-2.71 


-2.60 


-3.99 


-3.97 


-4.44 


-5.76 


-5.93 


-6.55 


-6.55 


-6.61 


-6.47 


-7.24 


-7.26 


-7.55 


0.73 


2.01 


0.74 


4.75 


0.05 


0.10 


6.52 


0.04 




0.0040 


2 


1.32 


-0.27 


-1.36 


-2.49 


-2.40 


-3.55 


-3.51 


-3.92 


-5.09 


-5.25 


-5.88 


-5.89 


-5.96 


-5.82 


-6.85 


-6.86 


-7.23 


0.91 


2.28 


0.76 


3.30 


0.08 


1.03 


3.12 


0.72 




0.0040 


3 


1.40 


-0.19 


-1.26 


-2.38 


-2.29 


-3.41 


-3.37 


-3.75 


-4.87 


-5.02 


-5.61 


-5.62 


-5.68 


-5.55 


-6.47 


-6.47 


-6.83 


1.00 


2.42 


0.77 


2.63 


0.10 


1.57 


0.42 


1.27 


|— — I 


0.0040 


5 


1.64 


0.06 


-1.04 


-2.23 


-2.13 


-3.27 


-3.21 


-3.58 


-4.66 


-4.80 


-5.35 


-5.36 


-5.40 


-5.31 


-6.09 


-6.09 


-6.38 


1.04 


2.47 


0.77 


2.39 


0.12 


1.85 


-0.93 


1.37 


w 


0.0040 


8 


1.67 


0.16 


-0.92 


-2.17 


-2.06 


-3.27 


-3.19 


-3.54 


-4.53 


-4.67 


-5.10 


-5.11 


-5.10 


-5.13 


-5.83 


-5.84 


-6.02 


1.11 


2.57 


0.77 


2.06 


0.14 


2.17 


-2.65 


1.64 


CO 


0.0040 


12 


1.76 


0.29 


-0.79 


-2.06 


-1.96 


-3.13 


-3.05 


-3.37 


-4.31 


-4.44 


-4.81 


-4.81 


-4.78 


-4.84 


-5.19 


-5.20 


-5.05 


1.15 


2.63 


0.77 


1.84 


0.15 


2.43 


-3.56 


1.80 


0.0040 


17 


1.72 


0.35 


-0.71 


-1.96 


-1.85 


-3.00 


-2.92 


-3.24 


-4.17 


-4.30 


-4.66 


-4.66 


-4.63 


-4.69 


-4.94 


-4.96 


-4.68 


1.18 


2.67 


0.77 


1.76 


0.17 


2.56 


-3.77 


1.80 




0.0080 


1 


1.49 


0.29 


-0.88 


-2.39 


-2.27 


-3.83 


-3.79 


-4.28 


-5.67 


-5.85 


-6.56 


-6.57 


-6.64 


-6.47 


-7.38 


-7.39 


-7.73 


0.77 


2.15 


0.81 


4.34 


0.07 


0.65 


5.78 


0.81 


C 


0.0080 


2 


1.76 


0.16 


-0.96 


-2.27 


-2.16 


-3.58 


-3.49 


-3.91 


-5.08 


-5.23 


-5.89 


-5.90 


-5.97 


-5.84 


-6.84 


-6.85 


-7.23 


0.99 


2.51 


0.82 


2.95 


0.12 


1.76 


1.33 


2.05 


0.0080 


3 


1.88 


0.27 


-0.84 


-2.15 


-2.05 


-3.46 


-3.37 


-3.78 


-4.89 


-5.04 


-5.64 


-5.65 


-5.70 


-5.61 


-6.54 


-6.55 


-6.91 


1.07 


2.64 


0.83 


2.44 


0.15 


2.24 


-1.24 


2.64 


0.0080 


5 


2.20 


0.69 


-0.41 


-1.81 


-1.69 


-3.22 


-3.11 


-3.53 


-4.63 


-4.78 


-5.35 


-5.36 


-5.39 


-5.34 


-6.17 


-6.18 


-6.51 


1.09 


2.65 


0.82 


2.27 


0.16 


2.49 


-2.51 


2.68 


r: 

:- 


0.0080 


8 


2.15 


0.75 


-0.31 


-1.76 


-1.65 


-3.38 


-3.24 


-3.69 


-4.74 


-4.90 


-5.40 


-5.42 


-5.40 


-5.46 


-6.03 


-6.06 


-6.31 


1.17 


2.82 


0.85 


1.94 


0.18 


2.85 


-4.06 


3.07 


0.0080 


12 


2.18 


0.87 


-0.16 


-1.63 


-1.52 


-3.30 


-3.16 


-3.61 


-4.64 


-4.80 


-5.26 


-5.28 


-5.24 


-5.33 


-5.75 


-5.78 


-5.96 


1.23 


2.92 


0.86 


1.70 


0.20 


3.13 


-5.19 


3.32 


X 


0.0080 


17 


2.23 


0.98 


-0.01 


-1.48 


-1.37 


-3.13 


-3.00 


-3.46 


-4.49 


-4.65 


-5.10 


-5.11 


-5.07 


-5.16 


-5.44 


-5.47 


-5.62 


1.27 


2.98 


0.86 


1.53 


0.21 


3.32 


-5.79 


3.41 


O 


0.0200 


1 


1.82 


0.69 


-0.36 


-1.89 


-1.77 


-3.61 


-3.53 


-4.05 


-5.44 


-5.62 


-6.41 


-6.42 


-6.50 


-6.32 


-7.36 


-7.36 


-7.74 


0.83 


2.34 


0.85 


3.83 


0.12 


1.53 


4.25 


2.80 




0.0200 


2 


2.31 


0.80 


-0.26 


-1.65 


-1.56 


-3.61 


-3.43 


-3.97 


-5.07 


-5.24 


-5.90 


-5.92 


-5.96 


-5.92 


-6.89 


-6.91 


-7.31 


1.04 


2.74 


0.87 


2.76 


0.17 


2.47 


-0.50 


4.01 


b 


0.0200 


3 


2.54 


1.00 


-0.04 


-1.50 


-1.41 


-3.72 


-3.50 


-4.08 


-5.08 


-5.26 


-5.79 


-5.83 


-5.82 


-5.91 


-6.61 


-6.69 


-7.02 


1.13 


2.94 


0.89 


2.29 


0.20 


3.02 


-3.33 


4.78 


1-3 


0.0200 


5 


2.76 


1.26 


0.24 


-1.20 


-1.11 


-3.49 


-3.26 


-3.84 


-4.84 


-5.02 


-5.55 


-5.59 


-5.57 


-5.68 


-6.24 


-6.33 


-6.65 


1.18 


3.00 


0.89 


2.01 


0.22 


3.35 


-4.67 


5.08 


o 


0.0200 


8 


2.57 


1.35 


0.39 


-1.01 


-0.94 


-3.38 


-3.14 


-3.74 


-4.74 


-4.92 


-5.44 


-5.48 


-5.46 


-5.58 


-6.05 


-6.16 


-6.46 


1.24 


3.13 


0.91 


1.74 


0.25 


3.69 


-5.82 


5.49 


p^ 


0.0200 


12 


2.47 


1.48 


0.57 


-0.81 


-0.75 


-3.35 


-3.09 


-3.71 


-4.70 


-4.90 


-5.42 


-5.47 


-5.45 


-5.60 


-5.95 


-6.11 


-6.40 


1.30 


3.26 


0.93 


1.51 


0.27 


4.00 


-6.80 


5.85 


c. 


0.0200 


17 


2.51 


1.57 


0.69 


-0.67 


-0.61 


-3.30 


-3.02 


-3.67 


-4.65 


-4.86 


-5.38 


-5.44 


-5.42 


-5.58 


-5.90 


-6.10 


-6.39 


1.34 


3.35 


0.94 


1.33 


0.29 


4.20 


-7.50 


6.12 


0.0500 


1 


2.22 


1.08 


0.17 


-1.16 


-1.09 


-4.45 


-4.16 


-5.03 


-6.21 


-6.43 


-7.18 


-7.24 


-7.27 


-7.33 


-7.81 


-7.96 


-8.24 


0.93 


2.94 


1.07 


3.25 


0.18 


2.45 


1.62 


5.76 


s. 


0.0500 


2 


2.80 


1.32 


0.32 


-1.04 


-0.97 


-4.14 


-3.86 


-4.72 


-5.77 


-6.01 


-6.65 


-6.75 


-6.77 


-6.95 


-7.35 


-7.60 


-7.92 


1.14 


3.23 


1.01 


2.42 


0.23 


3.32 


-3.70 


7.08 




0.0500 


3 


2.98 


1.48 


0.48 


-0.82 


-0.77 


-3.85 


-3.57 


-4.39 


-5.39 


-5.62 


-6.21 


-6.31 


-6.31 


-6.51 


-6.90 


-7.15 


-7.47 


1.22 


3.36 


0.99 


2.04 


0.27 


3.81 


-5.52 


7.89 




0.0500 


5 


2.78 


1.71 


0.75 


-0.57 


-0.52 


-3.66 


-3.38 


-4.22 


-5.21 


-5.47 


-6.02 


-6.14 


-6.14 


-6.40 


-6.62 


-6.96 


-7.29 


1.31 


3.52 


1.02 


1.71 


0.30 


4.32 


-6.97 


8.58 




0.0500 


8 


2.63 


1.87 


0.94 


-0.36 


-0.31 


-3.52 


-3.25 


-4.11 


-5.11 


-5.38 


-5.93 


-6.06 


-6.06 


-6.35 


-6.49 


-6.88 


-7.21 


1.35 


3.58 


1.02 


1.52 


0.32 


4.61 


-7.68 


8.95 




0.0500 


12 


2.42 


1.85 


0.98 


-0.26 


-0.22 


-3.35 


-3.08 


-3.92 


-4.92 


-5.18 


-5.73 


-5.86 


-5.86 


-6.13 


-6.28 


-6.65 


-6.99 


1.41 


3.69 


1.02 


1.32 


0.36 


4.96 


-8.52 


9.69 




0.0500 


17 


2.38 


1.91 


1.08 


-0.10 


-0.06 


-3.14 


-2.87 


-3.71 


-4.71 


-4.97 


-5.52 


-5.64 


-5.64 


-5.90 


-6.05 


-6.41 


-6.76 


1.45 


3.73 


1.01 


1.16 


0.38 


5.22 


-9.01 


10.28 





